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(57) A method of tracking a system (1 ) having sen- 
sors (3) (e.g. gas turbine engines) by providing a com- 
puter model (4) which, from changes in sensor readings, 
can track (follow) changes in performance parameters 
of the system. These may be e.g. efficiencies of various 
of components such as compressors and turbines. In 
the common embodiment, the method used (Singular 



Value Decomposition) is a way of computing the most 
likely solution when there are more performance param- 
eters than the number of sensors used and thus a de- 
gree of redundancy. The method allows system diag- 
nosis and most importantly a tracked model allows the 
selection of a tailored control regime (9) which results 
in optimum efficiency for the individual system lengine) 
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to tracking the performance of gas turbine enaines it .hliH nn . '^'^ specification refers mainly 
for any sys.em which incorporates mrsurerm e^^^^^^^^^ 

Such tracking models have practica^use TtooT. ,n 1^ exclusively applicable lo gas turbine engines 

control regimes can be app^d ole mo^T^ch ,h^^^^^^^^ °"T" performance of ,he engines, whereby varying 

determined. The op.imumTon.ro regTmrcan hen b ao^^^^^^^^^ - 

the state of Ihe engine in terms of t L2,lZr,^^ ^ '"^'"^ ^°^ds. depending on 

demand cnteria wiS var TZ^lTo'Ze^^^^^^^^ ^^^'"^^ P^^'^^' °P-aUng 

optimum control regime as peHtce chafes S^^ ^ ' ''""^ "^^^ "^^^ '° '^^ 

regime which achieves a number oftaranTeed oS^^^ f this standard engine is used to determine a control 

var Ji^^rtntThtdl^duTe^ngTnn:^^^^^^^^^^^ ^^^^'^ ''-'^ ^'"--es and tolerance 

measures such as specific fueJ consur^ tl o d dinrp^ Pe^'o^-^ance degrades causing performance 

causes each engine to changed aTr^Hala Lnse AnI?'"' ''''T'^" "^""^^ '^^-'-e 
soakage (rematching of the en Jne^^e 'o rrmaTchanoell b 17^ P^'^—- change ,s hea, 

and turbine efficiencies (these latter effeSs arl rSl^.^hE.^ f l '" '^'^^^^"^es) which affects compressor 

made when designing cLtroSerrfor a ^^0^!^^^^^^^^ Inevitably therefore, a number of compromises have .0 be 
the actual engine and the flx^modeUendt b^^^^^^^^ Z '"r""' ^"f^^^"^" 'between 

losses e g ,n fuel consumptSn ^en use^L,^^^^^^ "'""?"!^ '^^^ '° ^'9"'"^-"' 

-^-^ If optimum control is obtained for the indMduareng™^^ ' ""'"'^ corresponding gains can therefore be made 

.he elgi:roSrg^,;lrbTa,::;trut^^^^^^^ r "^-r^ ^ - -^.ne varia„on during 

that control optimisatL usina a va^n. moH . ^ °' """^'^'^ ^"9'"^ '=°"'^°l ^a'a. It ,s estimated 

foel consumpL . 1^ oI?of a oT^Trl ^rra^T'fbeS^^ ^"^'"^ ^ ^ ^ ^P-'^^ 

of this scale would be costly and hard won ihroMnh n?h ^""^ V °' temperature. Gains 

design. In fact additions 1, ol and hea^ ZlT '^^^'°P'^^'^'^' approaches such as improved turbor^achinery 

.ionassmai,aso.iorarnr:rror^^^^^^^^^^^ 

- for an aircraft when c uTgt a, iZTs k^^^^^^^^^^ ForeLmple' 
such systemsareSt I he 

Process for the F10O Engine" by J S O me S pub^^^^^^^^^ ^'''"^ 
" AIAA-92.3745. and ,n NASA te^hni.: r„ ^'^ V? T'""''' '"^'''"'^ °' Aeronautics and 

Astronautics 

Estimation Technigues Using Ka.man pTer'S:! ues '^^^ H liT m^erf te^' r"'°'" '"^'"^ 
performed on a computer model of the enainp and T ,1 ■ ^'^e optimisation described iherein ,s 

changes in sensor reading obtained fZ ih! . '"^'"^ '"''^ "'"^ °' '^^^^ '^^<=^"^9 "^^^-^^^ ,s ,0 use 

component performance 1 e to c^S ^^0?^,, h'h' ''"'""^ ^^'"'^ '° ^^'""^'^ - -gine 

« -after referredtoas S formanc^^^^^^^^ ^^'^^ ^'^ -"erna„voiv .nci here 

Of turbines or compressorchanaesr^-^^^ Parameters are for example efficiencies or flow c.pacties 

form Of correction terms which whTinputT^^^^^^^^^^^ ' ^^e .h 

changes as those readinos from ,hJ=^, ? '^^'^P^'^'ion o' identical modo^ sorter output 

said .0 match oMr ck the e gTe sucS T''"' '^^ ^^--^ ^---dJ: 

=- Which typically in addition to the rfn.T.^n , "'"^"^ ^ ''^'■""'^ thermodynamic mode- ... -c .name 
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sonable level of confidence that the model is a good representation of the engine This means that the sensors used 
must be widely distributed in terms of their independence 

Tracked models can be exploited practically to obtain control data which will maintain optimal efficiency for a certain 
required performance e g. specific fuel consumption in addition a knowledge of the change in performance parameters 
is useful in monitoring degradation of components and to investigate suitable maintenance action. 

The success of the model used for optimisation is crucially dependent on how well the model matches the engine 
since, as mentioned, the performance of every engine is different because of manufacturing tolerances, and will in any 
case deteriorate throughout the engine's operational life. Engine performance parameters are not directly measurable 
with engine control instrumentation but changes in their value can be estimated using prior knowledge of how such 
changes thereof affect changes in engine sensor outputs at a particular operating point. In mathematical terms, vari- 
ations in the set (or vector) of engine performance parameters which may occur (hereinafter denoted as dP)> gives 
rise to changes in the set (or vector) of engine sensor outputs (hereinafter denoted asdx),. when control demand inputs r 
to the engine are constant, i.e. at a particular operating condition. It ^s assumed that dP and dx are related linearly at 
a given engine operating point by a matrix of sensitivity coefficients C. such that:- 
dx = CdP 

dx contains m elements, m being the number of sensors used by the model, dP contains n elements, n being the 
number of performance parameters under consideration. C is a matrix of real numbers with m rows and n columns. 

NASA Technical Memorandum 104233 describes turbofan engine deterioration estimation by providing a tracking 
model of the engine which uses Kalman filtering techniques to determine dP from dx. There is a significant problem 
with this technique in that the tracking can only be performed if there are as many sensors used as performance 
parameters to be estimated i.e. if m = n. This forces various assumptions to be made, and results in the tracking not 
being exact. 

It is often the case that the number of performance parameters liable to measurable variation exceeds the number of 
sensors. A given set of changes in the sensor outputs of the engine at a particular operating condition could be ac- 
counted for by a variety of different performance parameter changes. However, although there may be no unique 
solution for the performance parameter variations causing a particular set of sensor variations, there is a component 
which is common to all possible solutions: the shortest component in vector space. The inventors have determined 
that this component gives the most likely and accurate solution. This is obtained by calculating a "pseudoinverse" of 
C which when obtained^ allows the most accurate dP estimate to be obtained, denoted dPe. 

It is an object of the present invention to provide a method of tracking the performance of a system, which allows 
performance parameters to be ascertained accurately, and further without the need for the number of sensors used 
by the tracking method equal to the number of performance parameters. 

The invention consists of a model of a system having input from one or more sensors of said system and having* 
means to relate the vector of changes in performance parameters of said system dP. to the vector of sensor output^ 
changes c|x. where the vectors are of different dimensions, and to determine that component of performance parameter 
changes which is common to all possible solutions and using this to determine performance parameters. 

The invention also comprises a method of tracking a system having sensors, comprising the steps of: 

a) providing a real time model having input from one or more of said sensors: 

bl storing at least one non-square coefficient matrix C relating sensor changes to performance parameter changes: 
and 

c) calculating a pseudoinverse of matrix C. where dx = C dP. where dx is the vector of system sensor changes 
and dp is the vector of performance parameter changes. 

b) from step (c) determining changes in performance parameters. 

The pseudoinverse of C is preferably calculated using a matrix method incorporating the technique of Singular 
Value Decomposition (SVD). which is a mathematical process used to analyse systems of linear simultaneous equa- 
tions, and will be described in more detail hereinafter. 

By way of example an embodiment of the invention will hereinafter be described with reference of the drawings of 
which' 

Figure 1 shows a block diagram of the tracking system for a gas turbine engine incorporating a tracking model. 
Figure 2 shows an expanded representation of figure 2 to show more detail of the tracking model 
Figure 3 shows graphs of the time response in determining performance parameters obtained in tests 
Figure 4 shows an optimising system incorporating a tracking model used to optimally control an engine 
Figure 1 shows a schematic figure of the tracking process used lo track the characteristics of a gas turbine engine 
( 1 i A thermodynamic model (4', is provided which is designed to track the aeroengine Input to the model are the values 
of the engine demand parameters from input demand parameter sensors (2) These represent the engine control inputs 
and engine operating environment parameters table l below shows a list of these parameters which are equivalent 



for both engine and model 
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Table 1 


Input demand parameters 


Variable 


Description 


wfdm 


Fuel supplied to the engine 


gvdm 


Engine inlet guide vane position 


sndm 


Engine nozzle actuator position 


hpbodm 


HP blow-off valve position 


Pti 


Total pressure at engine inlet 


ttl 


Total temperature at engine inlet 


pso 


Ambient pressure 



30 



25 



30 



35 



45 



SO 



uses easing engine conM ZITSZ.S!^ ' °' 

c»e, a. n,eL„ie„,s ^^:^z::::iZ'::z7L't,z:::':s:'^i::r ~ 



Table 2 



Output sensors 




Description 


Parameter 


nh 


high pressure spool speed 


nl 


high pressure spool speed 


PS 


low pressure compressor static pressure 


PS21 


by pass static pressure 


P21 


by pass total pressure 


P3 


compressor outlet pressure 


Ipsot 


low pressure stator outlet temperature 



parameters such as emc.encv and to^cLTrftv . h ''^"^'^ pararrielers Performar.ce 

...s enne .„g„e ,„ en, :;irr;:n-r3S==^^^^^^ 

Table 3 



Er^gine performance parameters (degradation) 



Variable 



^LPCI 
^LPCO 
®HPC 
^HPT 



^LPT 

^PCI 

^LPCO 

^HPC 

^HPT 

^PT 



Description 



increment in efficiency of Ip compressor (inner) 
Increment in efficiency of Ip compressor (outer) 
Increment in efficiency of hp compressor 
Increment in efficiency of hp turbine 
Increment in efficiency of Ip turbine 
Increment in flow capacity of Ip compressor (inner) 
Increment in flow capacity of Ip compressor (outer) 
Increment in flow capacity of hp compressor 
Increment n flow capacity of hp turbine 
Increment in flow capacity of Ip turbine 
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of changes m performance parameters on sensor outputs depend on the particular operating condition and have to be 
determined tn advance This will be explained toward the end of the specification 

The number of performance parameters n exceeds the number of sensors m. and equation 1 has no unique inverse. 
In other words, more than one set of performance parameter changes can be found which results in a particular set 
of output sensor changes. In order to determine the component which is common to all possible solutions and which 
gives the most accurate estimate, the pseudotnverse of C has to be found. As isshowrj^^y^iggre S.^liQiJ^ackgi^ cnocl^l 
stores previous sensor readings x tat a particular operating point) and ^CQm|)ute|^^e frbm previptjs^ engine 
output sensor readings stored by the model, x and the current engineToutptjr 

vecioYI^. This vector of sensor differences (or changes), dx undergoes three separate mathematical processes as 
shown. These are the Singular Value Decomposition (5). gain by factor K. (6) and signal integration with respect to 
time (7). Singular Value Decomposition (SVD) enables an approximation of a vector of engine performance variation 
elements. dP. to be obtained from a vector of engine output sensor changes output differences, dx. Singular Value 
Decomposition manages the redundancy in computing dP so as to calculate a "pseudoinverse" of matrix C. 

Singular Value Decomposition is a mathematical process which can be is applied to linear simultaneous equations 
represented by equation 1 . the exploitation of which enables accurate tracking to be achieved by accurately computing 
dP . The matrix C first undergoes decomposition. The Singular Value Decomposition for C is as follows: 
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25 



C =U S 0 ,v 

mxn mxm mxm mx(n-m) nxn 

Where U is a real mxm unitary matrix (representing a multidimensional axis transformation), the colunnns of which 
form an orthogonal basis for the m dimensional sensor space of C. Each of the columns is also scaled so that it has 
unit length and UU''' =: t^. V is an nxn unitary matrix (representing another multidimensional axis transformation) and 
the columns of the V matrix form an orthonormal basis for the n dimensional degradation vector and VV""" = 1^ and S^ 
is an mxm positive, real, diagonal matrix with elements arranged in monotonically descending order (representing a 
set of scaling factors), i.e. 
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SI 


0 


0 


0 


0 


S2 


0 


0 


0 


0 


S3 


0 


0 


0 


0 


S4 



Sm 



45 



where SI > S2 > S3 >Sr^., >Sn > 0 and ar^ the singular values The matrix S is defined as [S^ Ol'^ i e the S^ matrix 
packed with zeros in order to fill it to dimension of mxn i e. 
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defines the process of Singular V^lue DecoZc^^Sion and ^nwl^" ^^'"^ compensator shown in figure 1 . ms 

Further reference ,o SVO is "NumerlcaTredSrs .he^ a^^^^^^^^^ 

Teukolsky and W T VeUeriing, Car^bridge urersi.y Prel 1 99^^ ^ " ^ '^'-"-V S A 

tracCrrrhirrorn trerp':;; ~ reirr^^^-^^ - - - 

allows pertorr^ance parameters to be va ed 0' ^^^^^^^^ °' '^^ -9-- -h.ch 

readings are recorded. The simulation ,s run ala^sZ^n nl J T changes in output sensor 

poin.) and the steady state measurrmenUesponse T °" '^^^"^'"9 °" °P-^^""g 

duced one at a time on the simulation SSanaes dxT" f?"" " Perlormance parameters idP] are intro 
collection of this data is processed b ma.hema cal corr^o t^J^^^^^^ T"' ^he 

the division of sensor output vector chTngesTd ^^^^^^^^ computation of C obtained by performing 

In addition to the Sinoular vJJZT . ^ PSrf°^"^ance parameter vector changes idPJ, 
sains K ,one for ^^Cl^^^^^:::^ '^^::7TT " — '^^ 

the final value of dP more quickly The inlrlr and lif. , '""^ processing d^ so as .0 achieve 

The following-^escribes two'tes s rSemons'ra e II XT'"'"' "^^^ 
.r^vention. The apparatus is as shown Tn gu'eTe "rt ^3 tnsTea^^^ "' '""^""^ ^'^^^^^'"S '° '^^ 

simulation was used instead. Unlike a real enqfne ZTJllll ""^'"^ ^ thermodynamic 

ance parameters to be changed at m resuITin sensor ""f "°" °' '"P"' P^'^^^"^" 

model ,0 compute (or estimate) pertormanc pa am^^^ '""^^ °' "^'^^ ^^^^ '^V '^^^ 'racking 

metrackingmodeltobe^udged. ?hetherm2yZrsir^^ 

confusion with the tracking model will be referred hereinafter to as ,h J engfne ' '° 

Ixama^l Demonslraiinn nf tr....|nQ a( a O dP^.^n p^,.. 

Table 4 



Parameter 



®HPC 
®HPT 
®LPT 



Set elements of dP set m engme 



-0.01 
-0 01 
00 



F.n^l values of dP, calculated by tracking sy stem" 

-0.0147 
-0,0175 
-0 0031 
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Table 4 (continued) 



to 



Parameter 


Set elements of dP set in engine 


Final values of dP^ calculated by tracking system 


®HPT 


0.0 


0.0053 


©LPT 


0 


0.0033 


*LPCI 


-0.02 


-0.0195 


^LPCO 


-0.03 


0.00006 


^HPC 


-0.025 


-0.0135 


^HPT 


-0.02 


-0.0185 


*LPT 


0 


-0.0047 



The results of the tracking are also shown in figure 3. showing the time response of the tracking process. The step 
in engine degradation occurred 2 seconds after the start of recording. The computed values of pertormance parameters 
reach a steady state within 20 seconds. 

The results show that the elements of dPe . the estimates performance parameter changes, are respectable esti- 
mates of degradation and certainly sufficient to be indicative of damage or thermal non-equilibrium (e.g. eLp^ and 
®LPCo) Singular Value Decomposition therefore allows successful quantification of degradation components from a 
given set of sensors. 

Example 2 Demonstration of tracking at non C design point. 

A second test is hereinafter described to demonstrate a tracking of the moc3l at an engine operating point which 
is away from the operating point for which C was measured. Both engine and mode! were set up at an operating point 
test point, with the variable nl = 84.63%. In this case the U, S and VR matrices fixed in the tracking subroutine were 
computed from a C matrix which was computed at the different running condition^ that at which the test descnbed 
above was run before {i.e. nl =79.5%). The dP vector of performance parameter changes was input to the engine as 
before and the tracking model was left to run. The calculated changes in performance parameters in the last column 
of Table 5 show that the model tracks the engine outputs quite well, indicating that tracking performance is not signif- 
icantly degraded by moving the engine operating condition away from the condition at which C matrix was measured. 
It indicates an inherent robustness in the tracking process and it will help to minimise the data storage requirement by 
enabling an increase in the interval between C measurement points. 



Table 5 



35 


Parameter 


Values elements of dP set in engine 


Final values of dPe calculated by tracking system 




®HPC 


-0.01 


-0.0149 




®HPT 


-0.01 


-0.0188 




®LPT 


0.0 


-0.0032 


40 


®HPT 


0.0 


0.0054 




©LPT 


0.0 


0.0034 




^LPCI 


-0.02 


-0.0192 




^LPCO 


-0.03 


0.00005 




^HPC 


-0.025 


-0.0142 


45 


^HPT 


-002 


-0.0184 




^LPT 


0.0 


-0.0049 



The results in general show that linear analysis involving Singular Value Decomposition enables very good esti- 
mations of engine performance parameter changes to be calculated. Convergence of model and engtne can be en- 
hanced by combining integral feedback action with the Singular V^lue Decomposition process. There is no restriction 
to the number of degradation effects which may be included in the analysis. However, the accuracy of estimation of 
the individual degradation influences depends on the number and nature of the engine sensors used. Depending on 
the complexity of the system different or additional sensors may be used. 

As mentioned C data will vary depending on operating point In order to enhance the accuracy the C data may 
be stored by the tracking model lor a plurality ot different operating points. Depending on the operating point of the 
engine the C data stored for the closest operating point to that of the running engine can be selected and used by the 
tracking model Alternatively C data may be ascertained for a particular operating point by interpolation techniques. It 
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sensors all o( which provide information for the model to track the system Inout dL^n^ T descnplion of engine 

- rraMratrrhe^ 

b) minimisation of engine hot-end temperature whilst maintaining thrust at a set value or. 

c) maximisation of engine thrust whilst observing engine structural limitations 

Claims 



15 



20 



25 



45 



2. 



A model of a system as claimed in claim 1 and wherein the means to determine the chances in r^Prform^n.. 
parameters of said system comprises means to store at least one non-square o fficL ma'r, C reIaZ ne 
^rma^nce parameter changes to sensor changes and means for calculating a pseudo;:;lrTe ^i c ^her^ ' 

dx IS the vector of sensor changes and dP the vector of performance parameter changes. 

svTtem l°om "JT " ' °' ' '^^'^9 "^^"^ '° ^^'---e the operating pom, of said 

system from said sensors, and selecting C matrix data dependent on said operating point. 

^si^^rvreromp^^s^t,:^ ^^'^ "^^^-^ ^^'^-'^'^ - ^ - 

A model as claimed ,n claims 1 to 4 additionally having means to apply integration to signals of sensor changes 
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6. Monitoring means for a system having sensors comprising a model as claimed in claims 1 to 5. and means to 
indicate performance parameters 

7. Optimisation means for a system comprising a first model as claimed in claims l to 5 and a second model having 
input thereto system sensor outputs and performance parameters from said first model and comprising means to 
input to said second model varying inputs representing system input demand parameters, means to determine the 
operational pertormance of said second model depending on all inputs into said second model, means to select 
that set of input demand parameters which results in optimum operational performance of said second model, and 
means to apply said selected input demands parameters to the system. 

8. Optimisation means for a system as claimed in claim 7 wherein said operational performance is dependent on a 
set of demand criteria which can be varied. 

9. A method of tracking a system having sensors, comprising the steps of: 

a) providing a real time model having input from one or more of said sensors: 

b) storing at least one non-square coefficient matrix C relating sensor changes to perfornnance parameter 
changes: and 

c) calculating a pseudoinverse of matrix C, where dx - C dP. where dx is the vector of system sensor changes 
and dP is the vector of performance parameter changes. 

b) from step (c) determining changes in performance parameters. 

10. A method of tracking a system as claimed in claim 8 additionally including, between steps (a) and (b) of claim 9. 
■ the additional steps of: 

(a) determine system operating point from said sensors, and 

(b) selecting C matrix data dependent upon said operating point. 

11. A method of tracking a system as claimed in claims 9 or 10. where step (b) uses the method of Singular Value 
Decomposition. 

1 2. A method of tracking a system as claimed in claim 1 1 . additionally including the step of applying signal integration 
to the dx vector. 

1 3. A method of monitoring a system comprising obtaining performance parameters from the methods of claims 9 to 1 2. 

14. A method of optimisation of a system comprising sensors comprising the steps of: 

a) tracking a system as claimed in claims 9 to 12: 

b) providing a tracked model therefrom: 

c) inputting to said tracked model varying parameters representing system input demand parameters: 

d) determining the operational performance of said model: 

e) selecting that set of said input demand parameters which results in optimum operational performance of 
the model dependent on pre-selected criteria: and 

f) applying said selected input demand parameters to the system. 

15. A method optimisation of a system as claimed in claim 14 additionally including the step of selecting optimum 
operational performance depending on system demand criteria. 
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Fig.4. 
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(54) Means and method for system performance tracking 



(57) A method of tracking a system ( 1 ) having sen- 
sors (3) (e.g. gas turbine engines) by providing a com- 
puter model (4) which, from changes in sensor readings, 
can track (follow) changes in performance parameters 
of the system. These may be e.g. efficiencies of various 
of components such as compressors and turbines. In 
the common embodiment, the method used (Singular 



Value Decomposition) is a way of computing the most 
likely solution when there are more performance param- 
eters than the number of sensors used and thus a de- 
gree of redundancy.. The method allows system diag- 
nosis and most importantly a tracked model allows the 
selection of a tailored control regime (9) which results 
in optimum efficiency for the individual system (engine). 



Fig.4. 
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